e02bdf

e02bdf © Numerical Algorithms Group, 2002.

Purpose

E02BDF Evaluation of fitted cubic spline, definite integral

Synopsis

[defint,ifail] = e02bdf(lamda,c<,ifail>)

Description

 
 This routine computes the definite integral of the cubic spline 
 s(x) between the limits x=a and x=b, where a and b are 
 respectively the lower and upper limits of the range over which 
 s(x) is defined. It is assumed that s(x) is represented in terms 
                                                _      
 of its B-spline coefficients c , for i=1,2,...,n+3 and 
                               i                       
                                                       _        
 (augmented) ordered knot set (lambda) , for i=1,2,...,n+7, with 
                                      i                         
 (lambda) =a, for i = 1,2,3,4 and (lambda) =b, for 
         i                                i       
   _   _   _   _                  
 i=n+4,n+5,n+6,n+7, (see E02BAF), i.e.,
 
                               q      
                               --     
                         s(x)= >  c N (x).
                               --  i i
                               i=1    
 
        _    _                                                   
 Here q=n+3, n is the number of intervals of the spline and N (x) 
                                                             i   
 denotes the normalised B-spline of degree 3 (order 4) defined 
 upon the knots (lambda) ,(lambda)   ,...,(lambda)   .
                        i         i+1             i+4
 
 E02BDF can be used to determine the definite integrals of cubic 
 spline fits and interpolants produced by E02BAF.
 

Parameters

e02bdf

Required Input Arguments:

lamda (:)                             real
c (:)                                 real

Optional Input Arguments:                       <Default>

ifail                                 integer  -1

Output Arguments:

defint                                real
ifail                                 integer